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Abstract 

The supernova paradigm for the origin of galactic cosmic rays can be tested using multifrequency observations of 
both non-thermal and thermal emission from supernova remnants. The smoking gun of hadronic acceleration in these 
sources can, however, only be provided by the detection of a high energy neutrino signal. Here we apply the theory of 
non-linear particle acceleration at supernova shocks to the case of the supernova remnant RX J 17 13.7-3946, which is 
becoming the stereotypical example of a possible hadronic accelerator after the detection of high energy gamma rays 
by the HESS telescope. Our aim is twofold: on one hand we want to address the uncertainties in the discrimination 
between a hadronic and a leptonic interpretation of the gamma ray emission, mainly related to the possibility of a 
statistical uncertainty in the energy determination of the gamma ray photons in the TeV region. On the other we 
want to stress how a kin^ neutrino telescope would break the degeneracy and provide evidence for efficient cosmic ray 
acceleration in RX J1713. 7-3946. A 3cr evidence would require about two years of observation. 
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1. Introduction 

Supernova remnants (SNRs) are the most plausible 
sources of galactic cosmic rays. The recent detection 
of multi-TeV gamma radiation from several SNRs makes 
the case stronger, especially when coupled with multifre- 
quency observations of the same remnants. However, de- 
spite all this progress, the evidence that SNRs are indeed 
the main contributors to Galactic cosmic rays remains cir- 
cumstantial. A smoking gun evidence of efficient acceler- 
ation of cosmic rays in these sources can only come from 
the detection of high energy neutrinos, resulting from the 
decays of charged pions within the source. 

The acceleration of cosmic rays at supernova blast 
waves is well described by the non linear theory of diffu- 
sive shock acceleration (NLDSA) (see i28ll for a review). 
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This theory allows us to calculate the spectrum and spa- 
tial distribution of cosmic rays accelerated at a supernova 
shock taking into account the dynamical reaction of the 
accelerated particles on the shock and, in its most recent 
version l.6», J6i,il7il . also the generation of magnetic field 
through streaming instabiUty induced by the accelerated 
particles ITs , 12 ] and the dynamical reaction of the am- 
plified magnetic field on the plasma. Such a dramatic im- 
provement in the quality of our theoretical approach al- 
lows us to finally compare theory with observations and 
make testable predictions for future observations. 

One of the most important recent breakthroughs in es- 
tablishing the SNR paradigm for the origin of cosmic rays 
has been the detection of narrow filaments of non-thermal 
X-ray emission in the direction of supernova remnant rims 
ifiol Is, US 25 , i^l ■ The filaments are most commonly in- 
terpreted as the result of severe synchrotron energy losses 
of ultra relativistic electrons which are forced to radiate 
a large fraction of their energy (in the form of X-rays) in 
a narrow region downstream of the shock. The required 
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strengths of the magnetic field are of order ~ 100-500jL/G 
downstream of the shock. Such large levels of magneti- 
zation might be the manifestation of magnetic field am- 
plification in the upstream region, as due to the stream- 
ing of cosmic rays [35, 13, ^. Not only the detection of 
large magnetic fields signals for efficient cosmic ray ac- 
celeration, but in turn large magnetic fields are needed for 
increasing the maximum energy of accelerated particles 
to the knee region il4ll . namely for efficient acceleration. 
Efficient cosmic ray acceleration and magnetic field am- 
plification are two sides of the same coin. 

On the other hand, it has been proposed that the inter- 
pretation of the narrow X-ray rims may be flawed ||32[1 : 
they could in fact result from damping of magnetic field 
in the downstream region. The emission region would 
be narrow because the field disappears, not because par- 
ticles lose energy eff'ectively. The situation is currently 
subject of active debate: the damping would naively re- 
sult in the appearance of filaments not only in X-rays, but 
also in the radio emission, and at present there is no evi- 
dence for such a phenomenon. In addition, the absence of 
magnetic field amplification (or a mitigation of the eff'ect) 
would reduce the maximum energy achieved by the accel- 
erated particles, _unless the shock configuration is quasi- 
perpendicular |23ll. An additional possibility has been 
proposed in 12011 : the magnetic field could be amplified 
downstream (and not upstream) because of the develop- 
ment of corrugations on the shock surface which result 
in eddies in which the magnetic field winds up and gets 
amplified. In order for this scenario to lead to large maxi- 
mum momentum of the accelerated particles, the field di- 
rection must be very inclined with respect to the shock 
normal in the upstream region. 

A powerful diagnostic tool for particle acceleration in 
SNRs is represented by the gamma ray emission in the 
GeV-TeV energy region. The spectrum of the gamma ray 
emission, its extension to high energies and the shape of 
the cutoff all provide precious information on whether the 
radiation is of leptonic or hadronic origin. In ll30ll we dis- 
cussed in detail the application of NLDSA to the case of 
the supernova remnant RX J1713.7-3946. If the HESS 
data on this remnant are taken at face value, then it is 
problematic to fit them in the context of a leptonic model 
for at least two reasons: 1) the spectral cutoff' expected in 
leptonic scenarios leads to a gamma ray emission which 
falls short of the highest energy data points measured by 



the HESS telescope; 2) in any case, fitting both the TeV 
and X-ray data requires the presence of a very large back- 
ground of infrared photons in the remnant. 

The second point raised above should be considered as 
a circumstantial evidence against a leptonic model, but 
by itself it does not suffice to rule out this class of mod- 
els. In principle the first point is rather solid, but experi- 
mental effects might weaken its significance: a statistical 
(and systematic) uncertainty in the energy determination 
aff'ects the steep gamma ray spectrum by making it look 
smoother than it actually is so that convolving the theoret- 
ical prediction with an estimate of the uncertainty in the 
energy determination makes the leptonic spectrum look 
more similar to the hadronic case. We address this issue 
in a quantitative way here. 

The hadronic interpretation appears to be more sound, 
but there are shortcomings in this case as well: first, the 
thermal X-ray flux expected from RX J17 13.7-3946 is 
larger than the observed radiation of synchrotron origin. 
Second, the number of electrons needed to explain the 
observations is about 10""* of the number density of ac- 
celerated protons, at odds with the e/p ratio observed at 
Earth at energy ~ 1 - 10 GeV. 

Both these points are however rather weak at the 
present time. The first is based on the assumption that 
electrons and protons share the same temperature down- 
stream. This condition is hardly achievable and in fact one 
can easily argue that the temperature of electrons should 
be much smaller than that of protons. On the other hand. 
Coulomb scattering might be sufficient to raise the elec- 
tron temperature to a level large enough to excite emission 
lines of heavy elements li3Ll . The second point is equally 
weak in that electrons might be accelerated at different 
stages of the SNR. Moreover, recent data from PAMELA 
111] and ATIC US] suggest that a substantial contribution 
to the observed electron spectrum at Earth might come 
from sources other than SNRs 122113711 . 

The safest way of proving or rejecting acceleration of 
hadrons in RX J17 13.7-3946, as well as in other rem- 
nants, is to search for neutrinos produced in the decays of 
charged pions. In this paper we apply the NLDSA model 
developed by (jsl 0], and previously used to describe the 
multifrequency spectrum of RX J 17 13.7-3946 [30], in or- 
der to calculate the expected neutrino flux from this rem- 
nant. 

Previous attempts at estimating the neutrino flux from 
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SNRs have been typically based on phenomenological ap- 
proaches, building on the assumption of a power law ap- 
proximation for the gamma ray spectrum and simple scal- 
ing relations between the gamma ray and neutrino spectra. 

In 1 8] the authors estimated a flux of 40 neutrino in- 
duced muons (and antimuons) per krn^ per year from 
RX J1713. 7-3946, using a power law for the y-ray spec- 
trum fyiE) cc £ based on CANGAROO observations; 
the maximum neutrino energy was assumed to be ~ 10 
TeV. Neutrino oscillations, absorption and location of 
the detector were not taken into account. Such effects 
were included by lfl9ll . where fyiE) oc E^^-^ was used, 
the power law index being inferred from HESS data (no 
cut-off energy was assumed in this approach). The au- 
thors obtain a flux of neutrino induced muons N, 



10 km yr . Again based on HESS data, 114 111 used two 



fi+fi 



different parametrizations for the y-ray spectra, a power 
law plus an exponential cutoff at iiy^max - 12 TeV and 
a broken power law (with a knee at 6.7 TeV), leading to 



Nfj+p - 4.8 and 5.4 km"^ yr"' respectively. An attempt 
to extract the proton spectrum from the HESS data on the 
gamma ray flux (assumed to be of hadronic origin), and 
compute the neutrino flux based on that, was done in ||39ll . 

In the present paper we carry out the calculations by us- 
ing our model of NLDS A which provides a self-consistent 
description of the acceleration of cosmic rays in the rem- 
nant. This leads to a shape of the gamma ray spectrum 
which is not a simple power law, due to the non linear 
effects induced by the dynamical reaction of the accel- 
erated particles and the amplified magnetic field. At the 
same time we also obtain self-consistently the neutrino 
spectrum, and we use it to derive the expected number of 
events in a km^ neutrino telescope. 

The paper is organized as follows: in §|2]we summa- 
rize the technical aspects of NLDS A and its application to 
particle acceleration in supernova remnants. We also de- 
scribe the calculations of the non thermal radiation from 
RX J1713. 7-3946, with special attention for the gamma 
ray emission in both the hadronic and leptonic scenario. 
We discuss in detail the possibility to use present and fu- 
ture observations to discriminate between the two, once a 
statistical uncertainty in the energy determination of the 
photon events is taken into account. We show that a lep- 
tonic scenario convolved with a gaussian distribution of 
the photon energies with 6E/E ~ 30% or larger leads 
to the impossibility to distinguish the leptonic prediction 



from the hadronic one, at least using the HESS data on 
RX J1713. 7-3946. The implications for future gamma ray 
telescopes are also briefly discussed. On these premises 
it is very important to aim at the detection of the associ- 
ated neutrino signal, whose intensity is calculated in § [3] 
where we also estimate the neutrino induced muon signal 
in a kirP neutrino detector We conclude in § |4l 

2. NLDSA and the non thermal emission of 
RX J1713.7-3946 

2.1. Spectrum of accelerated protons and electrons 

In NLDSA theory, the overall shock structure and the 
outcome of the particle acceleration process are inextrica- 
bly linked. When acceleration is efficient, the pressure of 
accelerated particles affects the shock dynamics, leading 
to the formation of a precursor, namely a region where the 
fluid velocity progressively decreases while approaching 
the shock from far upstream. At the same time the stream- 
ing of accelerated particles is responsible for the instabil- 
ities that lead to magnetic field amplification. In turn, the 
fluid profile in the precursor and the amplified, turbulent 
magnetic field, with the scattering it provides (and possi- 
bly the induced energy losses), determine the efficiency of 
particle acceleration and the resulting spectrum, including 
its high energy cutoff. 

The shock structure and the accelerated particle spec- 

) — \ 

trum are computed as in |30]: the basic structure of the 
calculation is the same proposed in [5J and |i6J, but with a 
crucial new aspect taken into account, namely the dynam- 
ical reaction of the self-generated magnetic field, which 
is included following the treatment of 1 16] and 1 17]. This 
means that the conservation equations at the shock and in 
the precursor are modified so as to include the magnetic 
contribution. The compression factor at the subshock, 
Rsub, and the total compression factor, R,o,, are deeply af- 
fected by this change, resulting in a decrease of the com- 
pression ratio in the precursor, Rtot/Rsub, as soon as the 
amplified magnetic field contributes a pressure compara- 
ble to that of the thermal gas upstream. This smoothening 
of the precursor reflects in spectra of accelerated particles 
which are closer to power laws, though the concavity typ- 
ical of NLDSA remains visible ] 17]. 

The normalization of the proton spectrum is an out- 
put of our non linear calculation, once a recipe for injec- 
tion has been established. Following ifisll . particles are 
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injected immediately downstream of the subshock. The 
fraction of particles crossing the shock surface which are 
injected in the acceleration process, T/inj, can be written as 



(1) 



Here ^ ~ 2 - 4 is defined by the relation pinj = ^ 
where /?th,2 is the momentum of thermal particles down- 
stream. ^ parametrizes the poorly known microphysics of 
the injection process, but pth,2 is an output of the problem: 
as a result, the injection efficiency is affected by the dy- 
namical reaction exerted by the accelerated particles and 
by the amplified magnetic field. 

Finally, the maximum momentum of the accelerated 
particles is determined following lfl4ll for the computation 
of the acceleration time in the presence of a precursor 

For protons we use as a prescription the equality be- 
tween the acceleration time and the age of the SNR: 



tacc(Pp.mdx) — tsi' 



(2) 



For electrons, energy losses can be important. Their 
maximum momentum pe^mux is determined by equating 
the acceleration time with the minimum between the time 
for energy losses and the age of the remnant. The loss 
time of electrons over a cycle of shock crossing needs 
to be weighed by the residence times, tr, upstream and 
downstream, so that the condition for the maximum mo- 
mentum, in the loss dominated case, can be written as: 



tr,l(p) + tr,2(p) 



tr.2(P) 



(3) 



where t; denotes the loss time, and the indexes "1" and 
"2" refer to quantities measured upstream and down- 
stream respectively. The residence times in the context of 
the non linear theory of particle acceleration can be writ- 
ten expHcitely (from Eqs. (25) and (26) of jUl). Eq. © 
must be solved numerically for Pe,ms,x, contrary to the 
case of acceleration in the test particle regime. However 
an approximate analytical solution, valid when only syn- 
chrotron losses are important, was also proposed in 03011 . 

As to the electron spectrum at the shock, fefi(p), this is 
easy to calculate for p Pe,max- In fact, at a given mo- 
mentum p, the slope of the electron and proton spectrum 
is the same, if one assumes that both species experience 
the same diffusion coefficient. What is unconstrained a 



priori is the relative normalization of the two spectra, Kgp, 
which can only be obtained by fitting the observations. 
This is reasonable since electrons do not exert any appre- 
ciable dynamical reaction on the shock. 

The spectrum of electrons at energies around and above 
Pe,mix, namely the shape of the cutoff is harder to calculate 
in the context of non-linear theory. Since the spectra we 
find for electrons at p < pe^max are not far from being 
power laws with slope ~ 4, we adopt the modification 
factor calculated by [44] for strong shocks in test particle 
regime. The resulting electron spectrum at the shock, in 
the loss dominated case, is: 



fe.o(P) = Kep fpfiip) 



\+Q.523{plp,,^,,Y 



-P'Ipl.a 



(4) 

What is important to notice in this expression is that the 
cutoff is not a simple exponential, a fact which reflects 
in the shape of the synchrotron spectrum radiated by the 
electrons, making it different from what assumed by most 
of the previous work on the subject. On the other hand, 
if the maximum momentum is indeed determined by the 
age of the remnant, then the cutoff shape is expected to be 
exponential. 

2.2. Magnetic field amplification and compression 

The turbulent magnetic field close to the shock can be 
enhanced by several physical processes. However, here, 
as in 1 30], we focus on the amplification due to resonantly 
excited streaming instability induced by cosmic rays ac- 
celerated at the shock. Resonant streaming instability 13511 
is likely responsible for most of the magnetic field ampU- 
fication in SNRs after the beginning of the Sedov phase 
1 7], while the non-resonant mode of the same instability 
1 12] is more effective at earlier times. When the predic- 
tions of linear theory are extrapolated to the non-linear 
regime of field amplification (which one is forced to do 
for lack of a better treatment), the resulting field strengths 
are in agreement with the values inferred by identifying 
the thickness of the X-ray filaments with the synchrotron 
loss length of the highest energy electrons. The strength 
of the magnetic field at the position x upstream, 6B{x), in 
the absence of damping, can be estimated from the satu- 
ration condition, that, for modified shocks, reads ILL 



P,y{x) = U{X) 



-3/2 



1 - u{xY 



4 Max) 



(5) 
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where pw{x) = 6B{x)^/{^jtpoUq) is the magnetic pressure 
normahzed to the incoming momentum flux at upstream 
infinity, U{x) = u(x)/u() and Ma a - uq/va with va the 
Alfven velocity at upstream infinity, where only the back- 
ground magnetic field, Bq, assumed parallel to the shock 
normal, is present. Eq. (|5]l correctly describes the effect 
of compression in the shock precursor through the term 
U{xy^^^. For the upstream temperature that we adopt in 
RX J1713. 7-3946 (see below), damping in the upstream 
region is expected to be negligible. 

The magnetic field downstream of the subshock is fur- 
ther enhanced by compression, according to; 



Bi - RsubBu 



(6) 



where Bi is the magnetic field immediately upstream of 
the subshock and we have used the fact that we are dealing 
with Alfven waves, and hence turbulence perpendicular to 
the shock normal. 

2.3. Computation of the radiation fluxes and spectra 

The flux of non-thermal radiation at y-ray photon en- 
ergies is computed as in ^3^. In that work we consid- 
ered both possible scenarios for the origin of high energy 
photons in RX J1713. 7-3946, namely a leptonic origin, 
through inverse Compton scattering (ICS) of ambient low 
energy photons by accelerated electrons, or a hadronic 
origin, from decay following nuclear collisions of rel- 
ativistic protons. 

In the hadronic scenario, the y-ray flux from 
RX J1713. 7-3946 is accompanied by a flux of neutrinos 
coming from the decay of charged pions that are produced 
in nuclear collisions together with neutral pions. We com- 
pute the fluxes and spectra in both channels using the ap- 
proximated expressions for the cross sections as provided 
inRef. 1 24]. 

The y-ray and the neutrino fluxes, <1>^'(£',) (/ = y, v), 
produced by p-p collisions from a source located at dis- 
tance d from Earth, can be expressed as follows: 



(1)0(£,) = 



And^ 



J drn(r) J 



dEpfpir, Ep) 



do-i(E,„Ei) 



dEi 



(7) 

In this expression, the apex is used to indicate the neu- 
trino flux that would have been seen at the Earth in the 
absence of neutrino oscillations. For gamma rays the two 



fluxes clearly coincide. Here n is the gas (target) density 
in the SNR and fp{r, Ep) is the distribution function of ac- 
celerated protons at the location r at a given energy Ep. In 
general both quantities depend on the location in the shell, 
but here we assume them to be constant in the region en- 
closed between the contact discontinuity and the forward 
shock and vanishing outside this region. Finally, do-j/dEj 
is the inclusive differential cross section for the produc- 
tion of particles of type /. This quantity is usually ex- 
pressed through the total inelastic p-p cross section, o-jnei, 
and the dimensionless distribution function Fi(xi, Ep) for 
secondaries, as 



dcTiiEp, Ei) 

^inel (Ep) 



Fi(^Xi,Ep) 



(8) 



where jc, = Ei/Ep is the fraction of proton energy trans- 
ferred to the secondary particle. For the functions Fj, 
which effectively enclose all the details of the hadronic 
processes involved, we use the analytical approximation 
derived in ll24ll on the basis of numerical simulations of 
p-p collisions with the publicly available code SIBYLL. 
The analytical formulae provide a very good description 
of the flux and energy distribution of secondaries for en- 
ergies above 0.1 TeV. As far as photons are concerned, 
Fy also includes the contribution of 77 meson decay, in ad- 
dition to that of 71°, with an overall accuracy of order a 
few %. The estimate of neutrino fluxes is slightly less 
accurate, because Fy only includes the decay of charged 
pions, while neglecting the contribution from /T-mesons, 
and therefore leading to underestimate the neutrino flux 
by about 10 %. At energies lower than 100 GeV, and 
down to the rest energy of the 7r-meson, we use the ex- 
trapolated formulae provided again by ll24ll . that should 
be accurate within 10 %. 

The flux and spectrum of ICS photons is calculated by 
using the exact kernel for ICS, with the full Klein-Nishina 
(KN) cross section |33]. The main target photon field con- 
tributing to emission in the high energy y-ray band is pro- 
vided by dust-processed infrared photons, described by a 
blackbody spectrum with temperature ~ 20 K [34.1 . Fol- 
lowing the instance of the IR-i-Optical photon background 
in the interstellar medium (ISM) we assume that the ratio 
of the optical to infrared energy densities remains ~ 20, 
while the energy density of IR light, Wjr, is left as a free 
parameter (in the ISM, Wm ~ 0.05 eV cm"^). 
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In order to compute the spectra of accelerated particles 
and the resulting emission, we need estimates for a num- 
ber of environmental parameters relative to RX J1713.7- 
3946, namely its distance (d), age (tsNR), expansion ve- 
locity (mo), and the values of temperature (Tq), density 
(no), and magnetic field strength (Bo), in the surrounding 
medium. In addition to these, the only free parameters 
of the model are: ^, entering Eq. ([TJ, and Kep, defining 
the ratio between accelerated electrons and protons. In 
the ICS scenario an additional free parameter is the above 
mentioned Wjr. 

The uncertainties on the various parameters and how 
they affect the results of our calculations are thoroughly 
discussed in [30]. Here we only summarize the values 
that we have found to provide the best fit to the mul- 
tifrequency data in both scenarios. Our adopted values 
of distance, age, expansion velocity and temperature are: 
d - 1 kpc, tsm = 1600 yr (consistent with the historical 
Chinese record of a supernova explosion in AD 393 |43]), 
Mo = 4300 km/s, Tq = 10^ K. Other parameters are, in 
the hadronic scenario: hq - 0.12 cm"^, Bq = 2.6 yuG, 
^ - 3.8; in the leptonic scenario, instead: no = 0.01 cm"^, 
Bo = 1.5/iG,^ = 4.1. 

The resulting values of magnetic field strengths, accel- 
eration efliciency and maximum energy of the accelerated 
particles are quite diff'erent in the two cases: the hadronic 
scenario entails eflicient acceleration, with a fraction of 
accelerated particles corresponding to about 10 and a 
maximum proton energy exceeding lO''* eV. In the lep- 
tonic case the fraction of accelerated protons drops to less 
than 10"^, corresponding to an energy conversion effi- 
ciency of about 2 %. The magnetic field downstream is 
of order 20 fiG, to be compared with B2 ~ 100 fiG for the 
first scenario (in agreement with the value inferred from 
the synchrotron loss length interpretation of the rim thick- 
ness). The ratio between the number density of acceler- 
ated electrons and protons at the shock turns out to be 
Kep ~ 10 if acceleration is efficient. On the other hand, 
inefficient acceleration, and the lower value of the mag- 
netic field associated with the leptonic scenario, would 
favour Kgp ~ 10"^, in agreement with measurements of 
the diffiise galactic cosmic rays. However, in order to fit 
both the X-ray and the gamma-ray fluxes in the context 
of this purely leptonic scenario a local energy density of 
infrared radiation ~ 24 times larger than the galactic av- 
erage is required. 
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Figure 1 : Hadronic versus leptonic scenarios for tlie origin of TeV emis- 
sion from RX J1713. 7-3946. In botli cases, the flux is computed for 
tlie best-fit values of the parameters, as specified in the text. The thick 
solid line represents the spatially integrated spectral energy distribution 
of photons from n'^ decay in the hadronic scenario, while the thin solid 
line refers to ICS in the leptonic scenario. Symbols represent all avail- 
able HESS data and also plotted are EGRET upper limit and FermiLAT 
sensitivity for GeV energy photons from this source. Other Unes on the 
plot represent the result of convolution of the thin solid line with gaus- 
sians of different widths as explained in the text. 



In Fig. [U we plot the gamma ray flux as obtained 
through our calculations, for both a hadronic model (thick 
solid line) and a leptonic one (thick dashed line) (see 
1.30,1 ). Taken at face value, the curves in Fig. [T] clearly 
show that the hadronic scenario reproduces the high en- 
ergy observations much better, in particular the highest 
energy data points of HESS. On the other hand, the possi- 
bility to discriminate between the two models relies upon 
the assumption that the statistical uncertainty in the en- 
ergy determination is sufficiently small. In order to ad- 
dress this point, we calculate the theoretical prediction for 
the leptonic scenario in the case in which there is a statis- 
tical uncertainty in the energy determination cr(E) = xE 
with X = 0.2,0.3,0.4. Our results are shown again in 
Fig. [U (different line-types are labelled in the figure). 

The thin curves in the figure are obtained by convolving 
the predicted ICS flux at a given energy, (picsiE), with a 
gaussian of given width, namely : 



<D;cs(£)= dE' 



Oics(£0 
^J2n (tHE) 



exp 



{E-E'y 



'^2 



2o-2(£) 



(9) 



6 



The conclusion that these curves lead to is that a sta- 
tistical uncertainty larger than 30% in the energy region 
above 10 TeV would inhibit the discrimination between 
the two models, based on HESS data alone. The nominal 
energy resolution of HESS telescope is around 15% 
sufficient to allow for such discrimination. The main is- 
sue in this case is the data statistics: the main differences 
between the two scenarios, and the superior quality of the 
fit obtained within the hadronic one, show at the highest 
energy, where the statistical significance of the data points 
is not very high, as shown by the large error bars in Fig. [T] 
The role of the energy resolution is however important 
also in view of future gamma-ray telescopes, like CTA, 
whose design is currently being discussed by the scien- 
tific community. 

In addition to all this, it is worth keeping in mind that 
the Fermi satellite is expected soon to provide another 
crucial bit of information to this debate, in that the level of 
detection (or the non detection) should clarify the issue of 
a leptonic or hadronic origin for the gamma ray emission 
from RX J1713. 7-3946, independent of the shape of the 
high energy cutoff (see the low energy part of Fig. [T]). 



3. Neutrino signal 

The computation of the neutrino flux in the absence of 
neutrino oscillations was carried out in the assumption of 
perfect isospin symmetry (<l>;jO ^ (t>j^ ^ (t>n- ), which leads 
to equal neutrino and antineutrino fluxes of a given flavor 
For our purposes this is a good approximation. 

The neutrino flux at the Earth is related to through 
the oscillation probabilities: 



CD 



,(E,d)^ 2 P„,(E,d)<:>l,(E). 



(10) 



The transition probabilities, Pw, are in general func- 
tions of energy and travel length, but since at the ener- 
gies considered here the oscillation lengths are very short 
compared with the typical size of a SNR, the oscillation 
probability can be space averaged. The resulting flux of 

crossing the Earth, not taking into account absorption, 

is: 



with an identical equation holding for antineutrinos. The 
errors due to the uncertainties in the oscillation parame- 
ters are negligible (~ 5%). Since at the source and Vg 
are produced in a ratio {2 : 1 ), the effect of oscillations 
translates into a flux of muon neutrinos at Earth that is 
~ 50% of that produced at the source: Oy^ = 0.50^^. 

The flux of neutrinos and antineutrinos of each flavor 
expected at Earth is shown in Fig.|2] and compared with 
the flux of gamma rays (solid line). The flux of neutrinos 
has to be compared with the background whose main con- 
tribution comes from atmospheric neutrinos. The shaded 
region shown in Fig.|2]refers to the theoretical prediction 
for the atmospheric neutrino flux as we explain below. 

Following the estimates of 1 42 , 21 , 2l the atmospheric 
neutrino flux above 1 TeV can be approximated as: 



1 -2 -1 -1 

cm s sr 



0.4 <D° +0.2<^l 



(11) 



<Dj;;;,,(£y)^4.6xlO-«(^) "tcV 

(12) 

with an uncertainty of ~ 40% due to the experimental er- 
ror on the primary CR spectrum and composition and on 
theoretical models of hadronic interactions. Notice that 
the background flux also depends on zenith angle due to 
the different thickness of atmosphere to be crossed by cos- 
mic rays coming from different directions. This depen- 
dence is included in the thickness of the shaded region by 
averaging on all possible neutrinos' incoming directions 
(see [11] for a detailed description of the mean uncertain- 
ties in the neutrino flux determination). Notice that for 
atmospheric neutrinos we neglect the contribution of os- 
cillations, which is relevant only for energies ; 10 GeV 
L27.§6.1]. 

With respect to Eq. (fT2] i the atmospheric background 
shown in Fig. |2]is rescaled to a solid angle corresponding 
to a cone of semi-aperture 0.5 degrees. The motivation 
for choosing this value of the angle is twofold: it repre- 
sents a reasonable estimate for the angular resolution of 
a neutrino telescope at these energies, and it also corre- 
sponds to approximately the angular size of the shell of 
RX J1713. 7-3946. A smaller value of the aperture angle 
would imply a better chance to detect the signal than esti- 
mated below. 

In the following we specialize our predictions to the 
case of a kirv' neutrino telescope, for which the neutrino 
detection occurs by measuring the Cherenkov light from 
v-induced muons produced by charged-current interac- 
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Figure 2: Fluxes of photons and neutrinos (plus anti-neutrinos) from 
RX J1713. 7-3946 predicted in the hadronic scenario. Both electronic 
and muonic neutrinos contributions are shown. The shaded region rep- 
resents the range of theoretical predictions for the atmospheric muonic 
neutrino (plus anti-neutrino) flux integrated in a cone of semi-aperture 
0.5°. 



Figure 3: Neutrino event rate per year from RX J17 13.7-3946 (solid line) 
as a function of muon energy threshold for a neutrino telescope with 1 
km- elfective area. The dashed line shows the event rate expected from 
atmospheric neutrinos integrated in a cone of semi-aperture 0.5°. 



tions of Vfj and in the matter just below the detector. 
For the computation of the production rates of jj. and fi, 
we follow the method of Ref . il9tl . 

Since neutrinos can only be detected when the source 
is below the horizon of the detector, we introduce an aver- 
age live- time of the source, in the form of a parameter /uv, 
representing the fraction of time during which this condi- 
tion is satisfied. For a neutrino telescope as ANTARES 
and RX J17 13.7-3946 as the source, one can estimate 
/liv ~ 78% flof]. On the other hand, absorption of neutri- 
nos while crossing the Earth leads to an energy dependent 
reduction of the detected flux. In fact the Earth becomes 
opaque to neutrinos at energies Ey Z 1 PeV (when cr be- 
comes larger than lO""'-' cm^). For neutrino energies of 
~ 100 TeV (10 TeV) the signal is reduced by about 20% 
(5%). 

Taking into account these two effects, the number of 
muons with energy > £,/,, crossing an area A during 
the observation time T can be written as: 



= /livA T dEy (D,(£,) Y^(Ey, E„, 

Je,i, 



)x 



[1 -ay(Ey)] 



(13) 



with a similar equation holding foiNp. In Eq. ( fT3] ). ciyiEy) 
is the mean coefficient for neutrino absorption through 
the Earth, including only charged-current interactions and 



computed for a fixed Earth thickness resulting from av- 
eraging over the observation time lfl9ll . Yf^iEy, Eth) is the 
muon yield, namely the probability that a neutrino with 
energy Ey produces a muon with energy E^ > £,/, that 
crosses the detector area. This is obtained by integration 
over the muon energy, E^, of the neutrino interaction cross 
section multiplied by the muon range, R(Ef^, E,h): 

<A r ' dE^, "^{Ey, E^) R(E^, Erh) . (14) 
Je,„ dE, 

Here PhjO is the water density and Na is Avogadro's num- 
ber Both the muon range in water and the neutrino inter- 
action probability are taken from [19]. The latter is cal- 
culated by using the deep inelastic scattering formula for 
the charged-current cross-section, ctcc, and the distribu- 
tion functions of partons as calculated by [29]. 

In Fig.|3]the solid line shows the resulting A'^ + Np for 
A = 1 km^ and r = 1 yr as a function of the energy thresh- 
old. For E,h - 50 GeV the number of events per year is 
4.7, compatible with the findings of [41|. The shaded re- 
gion represents the muon background produced by atmo- 
spheric neutrinos, which can be computed as a function 
of the energy threshold, E,h, by substituting Eq. (fT2l i into 
Eq. ( fT3] l. The result is normalized to an aperture angle of 
0.5°, for the reasons explained above. The uncertainties 
represented by the shaded region have the same origin as 
for Fig. m 

The spectrum of neutrinos from RX J 17 13.7-3946 
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Table 1 : Compaiison between the expected signal and the atmospheric 
neutrino background. The last column shows the observation time (in 
years) required to obtain a signal from RX J1713. 7-3946 with a signifi- 
cance level of 3(T. 



£„,(GeV) 






yrs(3o-) 


100 


4.4 


4.7 


2.16 


500 


2.8 


1.3 


1.44 


1000 


1.9 


0.57 


1.41 



dominates the atmospheric background at energies above 
~ 300 GeV. In Table [T] we report the time (in years) of 
observation needed to obtain a signal with significance 
level of 3cr, assuming a simple Poisson distribution of 
events and a unit detection efficiency. One can see that if 
RX J1713. 7-3946 is indeed a hadronic accelerator a signal 
at the 3cr level could be seen in about 2 years of observa- 
tion at energies above 500 GeV. 

4. Conclusions 

In this paper we described the impact of non linear 
diffusive shock acceleration for the gamma ray and neu- 
trino production in SNRs. The intricate connection be- 
tween particle acceleration, magnetic field amplification, 
dynamical reaction of the particles and magnetic field, and 
the radiation produced by the accelerated particles (elec- 
trons and protons) in principle allow, if taken at face value, 
to impose strong constraints on the ability of SNRs to 
accelerate cosmic rays. An instance of how to use this 
powerful tool was provided in llsoll . where particle ac- 
celeration was described in terms of the non linear the- 
ory of Refs. [5i, ^]. The magnetic field amplification 
was assumed to be due to resonant streaming instabil- 
ity and both the dynamical reaction of the accelerated 
particles and of the amplified field were taken into ac- 
count. The basic conclusion reached in ll30ll is that the 
hadronic interpretation of the HESS data automatically 
leads to a downstream magnetic field which is in agree- 
ment with that inferred from X-ray observations if the 
rims of non-thermal emission are interpreted as the result 
of synchrotron losses. The non-thermal X-ray spectrum, 
as measured by Suzaku |38], was also reproduced with 
unprecedented accuracy. On the other hand, a satisfactory 



fit to the data within a leptonic scenario could only be 
achieved by considerably reducing the particle injection 
efficiency, through fine-tuning of the only free parameter 
of our calculations, ^, as defined in 32.11 Moreover, in 
order for the leptonic model to fit HESS data at all en- 
ergies one is forced to require the presence of a diffuse 
background of infrared light exceeding that observed in 
the interstellar medium by more than a factor 20. These 
requirements become less stringent if one decides not to 
include in the fit the highest energy HESS data points, 
which strongly constrain the maximum energy of the ra- 
diating electrons. 

Since the discrimination between the two models, 
hadronic and leptonic, is based on such a tricky region 
from the observational point of view, we decided to ad- 
dress here the issue of how well the photon energies need 
to be reconstructed at the telescope in order to tell the 
difference between the predictions of a leptonic and a 
hadronic model. The issue is of particular importance 
since the different predictions are of relevance in the en- 
ergy region where the particle spectra (and the gamma ray 
spectrum) are already sharply falling and a small statis- 
tical uncertainty in the energy determination may have 
sizable implications on the shape of the observed spec- 
trum. Our calculations show that a statistical uncertainty 
in the energy determination larger than about 30% would 
inhibit the possibility to discriminate between a hadronic 
and a leptonic interpretation of TeV data. Of course a sys- 
tematic uncertainty would strengthen the problem. These 
points need to be taken into account for the design of fu- 
ture telescopes, such as CTA. Needless to say that the ob- 
servation of TeV emitting SNRs with the Fermi gamma 
ray telescope will definitely contribute to settle the debate. 
This should certainly be the case for RX J1713. 7-3946 
(see also Q). 

In the absence of clear multifrequency evidence how- 
ever, the smoking gun that SNRs are efficient cosmic ray 
accelerators can only be provided by the unambiguous 
detection of neutrinos. Here we used the same non lin- 
ear theory of particle acceleration to infer the number of 
neutrino induced muons in a km^ neutrino telescope. The 
flux is compared with the appropriate background of at- 
mospheric neutrinos. At the distance of RX J1713. 7-3946 
and assuming a hadronic interpretation of HESS data, we 
predict that a 3cr detection should be achieved by hrr' 
neutrino telescope at the Antares' location (which is in 
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the right hemisphere to detect RX J1713. 7-3946) in about 
2 years of observation. 
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